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Abstract 

The first lecture provides an introduction to the physics of color 
superconductivity in cold dense quark matter. The main color su- 
perconducting phases are briefly described and their properties are 
listed. The second lecture covers recent developments in studies of 
color superconducting phases in neutral and /^-equilibrated matter. 
The properties of gapless color superconducting phases are discussed. 
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1 Introduction into color superconductivity 



1.1 Dense baryonic matter 

Almost all matter around us is made of several dozen chemical elements and 
their isotopes. Each element (atom) has a compact nucleus, made of protons 
and neutrons, and a cloud of electrons surrounding the nucleus. The proton 
and the neutron are baryons, i.e., "heavy" particles. They are about two 
thousand times heavier than the electron. (The actual values of the masses 
are m p « 938.3 MeV/c 2 , m n « 939.6 MeV/c 2 , and m c « 0.511 MeV/c 2 , see 
Ref. pp.) Thus, the mass of an atom comes mostly from its heavy nucleus. 
A typical size of a nucleus is about five orders of magnitude smaller than a 
typical atomic size. It is clear, therefore, that the density of matter inside a 
nucleus is much larger than the density of matter made of chemical elements. 
In fact, typical densities inside nuclei are of order 10 14 g/cm 3 . To get an 
impression how large this is, it suffices to say that a tablespoon of a material 
with this density would weigh about a billion tonnes. Creating even a tiny 
macroscopic sample of such matter in laboratory is impossible. On the other 
hand, baryonic matter at even higher densities exists in the Universe in the 
interior of compact stars. Therefore, it is of fundamental importance to study 
such very dense matter. 

Protons and neutrons, as well as all mesons and baryons that may appear 
in dense matter, are strongly interacting particles. Their properties and the 
properties of dense baryonic matter are, in principle, described by the micro- 
scopic theory of strong interactions, Quantum Chromodynamics (QCD). The 
only problem is that QCD is a rather complicated non-Abelian gauge the- 
ory. The Lagrangian density of the theory is written in terms of quark fields 
which carry color charges. Because of the property of confinement, however, 
quarks cannot exist in vacuum as free particles. Instead, they combine with 
other quarks and antiquarks to form color neutral hadrons (i.e., mesons and 
baryons) whose description is not easy in the framework of QCD. 

It is known that baryons are not point-like particles. They have a typical 
size of about 1 fm = 10~ 13 cm. In sufficiently dense matter, therefore, baryons 
can be forced to stay so close to one another that they would overlap. At 
such densities, constituent quarks are shared by neighboring baryons and, 
with increasing the density further, quarks should eventually become mobile 
over large distances. This means that quarks become deconfined. When this 
happens, it does not make sense to talk about hadronic matter any more. It 
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is transformed into quark matter. The properties of such matter is the topic 
of these lectures. 

It was suggested long time ago that quark matter may exist inside central 
regions of compact stars [2]. By making use of the property of asymptotic 
freedom in QCD [3], it was argued that quarks interact weakly, and that 
realistic calculations taking full account of strong interactions are possible 
for sufficiently dense matter The argument of Ref. jlj consisted of the 
two main points: (i) the long-range QCD interactions are screened in dense 
medium causing no infrared problems, and (ii) at short distances, the in- 
teraction is weak enough to allow the use of the perturbation theory. As 
we shall see below, the real situation in dense quark matter is slightly more 
subtle. 

1.1.1 Weakly interacting quark matter 

By assuming that very dense matter is indeed made of weakly interacting 
quarks, one could try to understand the thermodynamic properties of the 
corresponding ground state by first completely neglecting the interaction 
between quarks. In order to construct the ground state, it is important 
to keep in mind that quarks are fermions, i.e., particles with a half- integer 
spin, s — 1/2. Therefore, quarks should obey the Pauli exclusion principle 
which says that no two identical quarks can simultaneously occupy the same 
quantum state. 

In the ground state of zero temperature non-interacting quark matter, 
quarks occupy all available quantum states with the lowest possible energies. 
This is described formally by the following quark distribution function: 

f F {k) = e{ f i-E k ), at T = 0, (1) 

where fx is the quark chemical potential, and E k = y/k 2 + m 2 is the energy 
of a free quark (with mass m) in the quantum state with the momentum 
k (by definition, k = |k|). As one can see, /p(k) = 1 for the states with 
k < kp = V 'fi 2 — m 2 , indicating that all states with the momenta less than 
the Fermi momentum kp are occupied. The states with the momenta greater 
than the Fermi momentum kp are empty, i.e., /ir(k) = for k > kp. 

In the chiral limit, i.e., when the quark mass is vanishing, the pressure of 
zero temperature free quark matter is given by j3] 

P® = 2N f N c j —-- - E V ) Ofr-E k )-B= -^fi 4 - B, (2) 
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where the overall factor 2NfN c counts the total number of degenerate quark 
states, namely 2 spin states (i.e., s = ±1/2), Nf flavor states (e.g., up, down 
and strange), and N c color states (e.g., red, green and blue). The extra 
term B in the expression for the pressure (J2J), called the bag constant, was 
added by hand. This term effectively assigns a nonzero contribution to the 
vacuum pressure and, in this way, provides the simplest modelling of the 
quark confinement in QCD [UJ. [Note that the vacuum pressure is higher 
than the pressure of quark matter when fi < (4ir 2 B /Nf) 1 ^ 4 .] Here, we ignore 
the requirements of the charge neutrality and the (3 equilibrium in quark 
matter. These will be addressed in detail in the second lecture. The energy 
density of quark matter reads 

Thus, one arrives at the following equation of state: 

p(°) = I( e (°)-45). (4) 

This equation of state could be easily generalized to the case of nonzero 
quark masses and/or to the case of non-equal chemical potentials for differ- 
ent quarks. Also, it could be further improved by adding the lowest order 
corrections due to the interaction. For example, to leading order in coupling, 
the interaction results in the following correction to the pressure j7]: 

^ = - ^"V + 0(^), (5) 

where a s = g 2 /4:7t is the value of the running coupling constant of strong 
interactions defined at the scale of the quark chemical potential. The next 
to leading order corrections can also be calculated [Zj. 

It is generally believed that the equation of state in Eq. (J2J), or its gen- 
eralization, provides a good approximation for the description of weakly in- 
teracting quark matter. Below we argue, however, that some underlying 
assumptions regarding the ground state of such quark matter are not cor- 
rect. 



1.2 Cooper instability and color superconductivity 

It appears that the perturbative ground state of quark matter, characterized 
by the distribution function in Eq. is unstable when there is an attractive 
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(even arbitrarily weak in magnitude !) interaction between quarks. This is 
because of the famous Cooper instability [S] that develops as a result of the 
formation of Cooper pairs (g^ g_k) made of quarks from around the highly 
degenerate Fermi surface, i.e., quarks with the absolute value of momenta 
k ~ kp. Cooper pairs are bosons, and they occupy the same lowest energy 
quantum state at zero temperature, producing a Bose condensate. In the 
presence of such a condensate of Cooper pairs, the ground state of quark 
matter becomes a (color) superconductor. This is very similar to the ground 
state of the electron system in the Bardeen-Cooper-Schrieffer theory of low- 
temperature superconductivity j^J. The only real difference comes from the 
fact that quarks, unlike electrons, come in various flavors (e.g., up, down 
and strange) and carry non-Abelian color charges. To emphasize this dif- 
ference, superconductivity in quark matter is called color superconductivity. 
For recent reviews on color superconductivity, see Ref. [TU] . 

It was known for a long time that dense quark matter should be a color 
superconductor JH^|. In many studies, however, this fact was commonly 
ignored. Only recently, the potential importance of this phenomenon was ap- 
preciated. To large extent, this was triggered by the observation in Ref. J3| 
that the value of the color superconducting gap could be as large as 100 MeV 
at baryon densities existing in the central regions of compact stars, i.e., 
at densities which are a few times larger than the normal nuclear density, 
n ~ 0.15 fm~ 3 . This very natural estimate for the value of the gap in QCD, 
in which a typical energy scale itself is 200 MeV, opened a wide range of new 
theoretical possibilities, and the subject bursted with numerous studies. The 
main reason is that the presence of such a large energy gap in the quark spec- 
trum may allow to extract clear signatures of color superconducting states 
of matter in observational data from compact stars. 

As in low-temperature superconductors, one of the main consequences of 
color superconductivity in dense quark matter is the appearance of a nonzero 
energy gap in the one-particle spectrum, 



where A is the gap. The presence of the gap in the energy spectrum should 
affect transport properties (e.g., conductivities and viscosities) of quark mat- 
ter. Thus, if quark matter exists in the interior of compact stars, this will 
be reflected, for example, in the cooling rates and in the rotational slowing 
down of such stars. Also, a nonzero gap modifies thermodynamic properties, 




(6) 
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e.g., the specific heat and the equation of state. In application to stars, this 
could modify theoretical predictions for the mass-radius relations, or even 
suggest the existence of a new family of compact stars. 

Color superconductivity may affect directly as well as indirectly many 
other observed properties of stars. In some cases, for example, superconduc- 
tivity may be accompanied by the baryon superfluidity and/or the electro- 
magnetic Meissner effect. If matter is superfluid, rotational vortices would 
be formed in the stellar core, and they would carry a portion of the angular 
momentum of the star. Because of the Meissner effect, the star interior could 
become threaded with magnetic flux tubes. In either case, the star evolution 
may be affected. 

In general, it is of great phenomenological interest to perform a system- 
atic study of all possible effects of color superconductivity in compact stars. 
Before this can be done, however, one needs to know the structure of the 
QCD phase diagram and properties of various color superconducting phases 
in detail. Despite the recent progress in the field, such knowledge still remains 
patchy. While many different phases of quark matter have been proposed, 
there is no certainty that all possibilities have already been exhausted. This 
is especially so when additional requirements of charge neutrality and (5 equi- 
librium are imposed. 

In the rest of this lecture, the main color superconducting phases, that 
involve different number of pairing quark flavors, will be reviewed briefly. The 
role of charge neutrality and j3 equilibrium in quark matter will be discussed 
in the second lecture. 

1.3 Two-flavor color superconductivity (Nf = 2) 

It is instructive to start this lecture from the description of the simplest color 
superconducting phase, namely the two-flavor color superconductor (2SC). 
This is a color superconducting phase with the spin-0 Cooper pairing in quark 
matter made of up and down quarks. 

As we discussed in the previous subsection, color superconductivity comes 
as a result of the Cooper instability driven by an attractive interaction be- 
tween quarks. To study this instability in detail, one needs to specify the 
origin and the strength of the interaction. At very large densities, when 
quark matter is weakly interacting, the dominant interaction between quarks 
is provided by the one-gluon exchange. Then, the Cooper pairing dynamics, 
including the screening effects due to dense medium as well as the effects due 
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to higher order corrections, can be studied from first principles within the 
framework of QCD. Some key details of the analysis in this regime are given 
below. 

We consider quark fields in the fundamental representation of the 
SU(3) C color gauge group. These fields carry flavor (i = 1,2) and color 
(a = 1,2,3) indices. The QCD Lagrangian density reads 

£qcd = ft (vfdp + T V - mf } ) ft + gA A ft^T A ft - \g a u g a ^, (7) 

where A A is the vector gauge field in the adjoint representation of SU(3) C , 
G A V = d^A^ — d v A A + gf ABC A^A^ is the field strength, and the generators 

of color transformations are defined as = | (■^ j4 ) b where X A are the Gell- 
Mann matrices. The current quark masses and the quark chemical potential 
are denoted by mf and /i, respectively. At high densities, one can neglect 
small current masses of quarks. Then, the QCD Lagrangian density (|7|) 
becomes invariant under SU(2)ixSU(2)^ global chiral transformations. 

In QCD, the constituent masses of quarks rrii are generated dynamically, 
and they can be very different from the current masses mf\ appearing in 
the Lagrangian density. At zero density, for example, typical values of the 
constituent quark masses are of order |m n 313 MeV even in the chiral 
limit. At high densities, on the other hand, the masses of the up and down 
quarks become small. This is because they are proportional to the value 
of the chiral condensate (iPl^r) which melts in dense matter. Thus, it is 
often also justified to neglect the constituent quark masses in studies of color 
superconducting phases. 

The interaction vertex in the Lagrangian density (J2J) has a non-trivial 
color structure given by the color generators T^. As a result, the quark- 
quark scattering amplitude in the one-gluon exchange approximation is pro- 
portional to the following color tensor: 

^ T aa iT b , h = — — — (S aa i5b'b — bab'Sa'b) H TT7 — {Saa'&b'b + ^ab'^a'b) ■ (8) 

A=l 4iV c 4iV c 

The first antisymmetric term corresponds to the attractive antitriplet chan- 
nel, while the second symmetric term corresponds to the repulsive sextet 
channel, see Fig. ^ It is the first antisymmetric antitriplet channel that 
plays the crucial role in Cooper pairing. 
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= ( \ + 6/ 

-p -Jc Attractive 

Figure 1: The diagrammatic representation of the one-gluon exchange inter- 
action between two quarks in QCD. The color structure of the correspond- 
ing amplitude contains an antisymmetric antitriplet and a symmetric sextet 
channel. 




Here, it is appropriate to note that quark matter is unlikely to be truly 
weakly interacting at moderate densities existing in the central regions of 
compact stars, i.e., n<10n . In this case, the use of the microscopic theory 
of strong interactions is very limited, and one has to rely on various effective 
models of QCD. A very simple type of such a model is the Nambu-Jona- 
Lasinio (NJL) model with a local four-fermion interaction ^3]. Originally, 
this model was introduced for the description of nucleons with dynamically 
generated masses. Nowadays, it is also commonly used for the description of 
quarks. For example, the NJL type models were used in Refs. [TSI E3 EE Ej 
for the description of color superconducting phases. One of the simplest NJL 
models, that respects the SU(2)iXSU(2)^ global chiral symmetry (in the 
limit m| — > 0), is defined by the following Lagrangian density [T5] : 

+ G D (i^ c ee a ^)(i^ee a ^ c ), (9) 

where ip c = Cip T is the charge-conjugate spinor and C = «7 2 7° is the charge 
conjugation matrix. The matrix C is defined so that C^^C^ 1 = — 7T. Re- 
garding the other notation, f = (t^t^t 3 ) are the Pauli matrices in the 
flavor space, while (e) lk = e lk and [t a ) bc = e abc are the antisymmetric ten- 
sors in the flavor and in the color spaces, respectively. The dimensionful 
coupling constant Gs = 5.01 GeV -2 and the momentum integration cutoff 
parameter A = 0.65 GeV (which appears only in loop calculations) are ad- 
justed so that the values of the pion decay constant and the value of the 
chiral condensate take their standard values in vacuum QCD: F n = 93 MeV 
and (uu) = (dd) = (—250 MeV) 3 |18j . The strength of the coupling constant 



£njl = C(*7^ + 7 V-™H< + Gs 
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Go is taken to be proportional to the value of G$ as follows: Gp = r/Gs 
where rj is a dimensionless parameter of order 1. Note that the value of 77 
is positive which corresponds to having the antisymmetric diquark channel 
attractive. This is motivated by the microscopic QCD interaction, as well as 
by the instanton motivated models [T3] . 

1.3.1 Color and flavor structure of the condensate 

Now let us discuss the color and flavor structure of the condensate of Cooper 
pairs in the 2SC ground state. This can be determined from the nature of 
the diquark interaction and the Pauli exclusion principle. 

The color structure of the one-gluon exchange interaction in Eq. (jHJ) shows 
that the color antisymmetric antitriplet channel is attractive. The antisym- 
metric color channel is also attractive in the instanton motivated models of 
Ref. [T3j. This is reflected in the color structure of the condensate of spin-0 
Cooper pairs, 



which is antisymmetric in the color indices of the constituent quarks. As one 
can check (not shown explicitly), this spin-0 condensate is also antisymmetric 
in Dirac indices. Then, in accordance with the Pauli exclusion principle, it 
should be antisymmetric in flavor indices as well. The presence of the 7 s 
matrix on the left hand side of Eq. (|rUjl makes the condensate a scalar. A 
similar condensate without 7 s is a pseudoscalar which is disfavored in QCD 
because of the instanton effects. 

The color antitriplet condensate (jlUJl may have an arbitrary orientation in 
the color space. By making use of the global color transformations, this ori- 
entation can be changed as convenient. It is conventional to choose the con- 
densate to point in the third (blue) color direction, ((V' C '). 7 5 ^) ~ £ij£ ab3 - 
In this case, the Cooper pairs in the 2SC phase are made of the red and green 
quarks only, and the blue quarks do not participate in the pairing. These 
unpaired blue quarks give rise to gapless quasiparticles in the low-energy 
effective theory. 

The flavor antisymmetric structure in Eq. (|10|) corresponds to a singlet 
representation of the global SU(2) L xSU(2) fi chiral group. This means that 
the chiral symmetry is not broken in the 2SC ground state. In fact, there are 
no other global continuous symmetries which are broken in the 2SC phase. 
There exist, however, several approximate symmetries which are broken. One 




(10) 
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of them is the approximate U(l)^ symmetry which is a good symmetry at 
high density when the instantons are screened ^H]- Its breaking in the 2SC 
phase results in a pseudo-Nambu-Goldstone boson [20]. Additional four ex- 
otic pseudo-Nambu-Goldstone states may appear as a result of a less obvious 
approximate axial color symmetry discussed in Ref. [2*T] . 

In the 2SC ground state, the vector-like SU(3) C color gauge group of QCD 
is broken down to the SU(2) C subgroup. Therefore, five out of total eight 
gluons of SU(3) C become massive due to the Anderson- Higgs mechanism. 
The other three gluons, which correspond to the unbroken SU(2) C , do not 
interact with the gap less blue quasiparticles. They give rise to a pure SU(2) C 
gluodynamics. The red and green quasiparticles decouple from this low- 
energy SU(2) C gluodynamics because they are gapped |22j. 



1.3.2 Quark propagator in the Nambu-Gorkov basis 

To simplify the study of diquark condensates in dense quark matter, instead 
of using the Dirac spinors for quark fields, it is more convenient to introduce 
the following eight-component Nambu-Gorkov spinors: 



* = 



(ii) 



where ip is the charge-conjugate spinor, defined earlier. The structure of 
the inverse quark propagator in this basis reads 



o-l 

d 



where 



Go 







Gr 



-1 



7^ ± /ry 



121 



(13) 



are the inverse Dirac propagators for massless quarks (Gq) and charge- 
conjugate quarks (Gq ). 

In the ground state with Cooper pairing, see Eq. (JTUJ), the quark propa- 
gator in the Nambu-Gorkov basis also has nonzero off-diagonal components, 



s~ 



G r 



A- 



Gr 



-1 



(14) 
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where A = — ie'e^/S. and A + = 7 (A _ )^7° = —ie 3 £ / y 5 A* are the matrices 
in the Dirac space, and A is the diquark gap parameter. This is the simplest 
ansatz for the quark propagator in the color superconducting ground state. 
It can be made more general, for example, by adding a regular (diagonal) 
part of the quark self-energy. Such a self-energy plays an essential role in the 
dynamics of chiral symmetry breaking. In dense quark matter, on the other 
hand, it is not so important. Thus, for simplicity, it is neglected here (the 
effects of the diagonal part of the self-energy were studied in Ref. [2HD- 

From the inverse propagator in Eq. (|14|). one can derive the following 
expression for the quark propagator: 



S 




(15) 



where 



G ± 



(<*)"' 



-G^A ± G ± . 



A^GJA ± 



1 -1 



(16) 
(17) 



This structure of the quark propagator carries basic information about the 
properties of the 2SC ground state. For instance, the location of its poles in 
the energy plane gives the quasiparticle dispersion relations. 

The only quantity that remains unknown in the ansatz (|14|) is the value 
of the gap parameter A. This should be determined from a self-consistent 
gap equation that takes into account the pairing dynamics between quarks. 
This is discussed in more detail in the next three subsections. 



1.3.3 One-gluon exchange interaction 

In dense QCD, in which quarks are weakly interacting, the one-gluon ex- 
change is the dominant interaction between quarks. It appears, however, 
that the one-gluon interaction is partially screened by surrounding dense 
medium. Thus, before studying the pairing dynamics between quarks, we 
first address the modification of gluon properties due to dense medium. 

The inverse propagator of the medium modified gluon can be written in 
the following form: 

M:MM;>n- us) 
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where IT^f is the gluon self-energy, or the polarization tensor. The gluons 
with soft momenta, p^,g s ^, play the key role in Cooper pairing of quarks. 
The main contribution to the corresponding soft gluon polarization tensor 
comes from the quark loop with hard internal momenta of order \i. This 
contribution is large compared to all other contributions in QCD (e.g., those 
from the ghost and the gluon loops) because it is proportional to the density 
of quark states at the Fermi surface (as well as to the coupling constant 
a s ), i.e., Il^f ~ a s /i 2 . The corresponding approximation, in which the sub- 
leading terms suppressed by powers of the quark chemical potential /i are 
neglected, is called the hard dense loop (HDL) approximation. 

The calculation of the polarization tensor in the HDL approximation was 
performed in Refs. |2"H 125] (in the case of dense QCD in 2 + 1 dimensions, 
similar calculations were performed in Ref. [2E])- Here, we present the final 
result and briefly comment on its main features. 

The HDL polarization tensor can be written as follows: n^f = 5^11^, 
where the explicit form of the components of are given by 

n°> ,p) = Mp ,p), (19) 
n 0i (p ,p) = po^iMpo,p), (20) 

n*(Pb,p) = f^-^n t (p ,p) + ^4 n '(^p), (2i) 
\ p j p p 



with 



nz(po,p) = m 2 D 



f^ln 


Po+P 


- 1 


\2p 


Po-P 





2p 



(22) 



1 P 2 

n t (po,p) = -m 2 D - —Ui(p ,p). (23) 

Here, we use the notation p 2 = p 2 , P 2 = p\ — p 2 , and m 2 D = 2Nfa s /i 2 /n. It is 
remarkable that this result coincides with the polarization tensor derived in 
the framework of the classical transport theory of dense Yang-Mills plasma 
[2*7j . One can check that the HDL polarization tensor is transverse, i.e., 

PTV(p ,p) = 0. (24) 

It is convenient to write the tensor Tl^ in terms of two transverse projection 
operators, 

U, u = -OWn t + (2U t - m 2 D ) . (25) 
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where 



Oj$ = g„-u,u u + ^, (26) 

Of} = V*-^-^, (27) 

and Up = (1, 0, 0, 0) is the time-like four-vector that defines the rest frame of 
dense quark medium. In order to have a complete set of projectors such that 
Y,i Ofy = g^, one has to include also the following longitudinal projector: 

Of} = (28) 

By making use of this set of three projectors, inverse gluon propagator in 
Eq. (jTSjl can be represented as follows: 



V -l) AB = i5 AB 



p 2 - n t ) o« + (p 2 + 2n, - ml) o% + ^V) 



(29) 

where A is the gauge fixing parameter. This expression can be easily inverted, 
leading to the following result for the gluon propagator: 

K B - (j^V« + p2 + 2 n,- m ^ g + ■ < 3 °> 

The gluon modes which are transverse and longitudinal with respect to the 
three-momentum p are called the magnetic and electric modes, respectively. 
These two types of physical gluon modes are described by the following fields: 



a 



(mag) = Q m A ^ (31) 

< } = OftA". (32) 

From Eq. (J3U|) one can see that the propagators of these modes are naturally 
separated from each other, as well as from the unphysical longitudinal mode, 
described by the field a{p = 0$A V . 

In the study of the Cooper pairing dynamics, the most relevant part 
of the one-gluon exchange interaction is the interaction with the space-like 
momenta such that p <C p. This should be clear from the kinematics of 
the quark scattering around the Fermi surface: the energy exchange between 
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pairing quarks is typically much smaller than the change of their momenta. 
It is justified, therefore, to use the approximate form of U t , obtained in the 
region p <C p, 



By substituting this asymptote into the gluon propagator (J5UJ) . one can check 
that the magnetic modes produce a long-range (dynamically screened) inter- 
action, while the electric mode works only at short distances, r<l/m£,. 

In the next subsection, the HDL gluon propagator is used in the study of 
the Cooper pairing dynamics. As we shall see there, the long-range interac- 
tion mediated by the magnetic modes will play a particularly important role 
in the dynamics. 

1.3.4 Gap equation in QCD 

Color superconductivity is an essentially non-perturbative phenomenon which 
cannot be addressed with perturbative techniques even if the theory is weakly 
interacting. The standard method for studying superconductivity is the 
method of the Schwinger-Dyson equations. In general, the complete set of 
Schwinger-Dyson equations contains an infinite number of coupled equations 
for Green's functions. There exist no systematic ways of solving such equa- 
tions exactly. However, there exist various standard approximations that 
can be justified if the dynamics is controlled by a small parameter (e.g., a 
weak coupling constant, or a small 1/N parameter where N is the number 
of flavors or colors in the model). 

In the case of QCD at asymptotic densities, the main features of its 
dynamics are well described by the so-called improved rainbow approxima- 
tion. In this approximation, one uses the bare quark-gluon vertices in the 
Schwinger-Dyson equation for the quark propagator. One also uses the gluon 
propagator with the screening effects taken into account in the HDL approx- 
imation. (Note that, in the simple rainbow approximation, the screening 
is usually neglected.) The graphical representation of the Schwinger-Dyson 
equation is shown in Fig. |21 The analytical form of this equation reads 



S-\K) = S \K) + Ana s J j^Tp{P)Y B u V^ B {K - P), (34) 




(33) 
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Figure 2: The graphical representation of the Schwinger-Dyson equation in 
the improved rainbow approximation. The quark propagator in the Nambu- 
Gorkov basis and the gluon propagator are denoted by the solid lines and 
the wavy lines, respectively. 



where the quark-gluon vertices in the Nambu-Gorkov basis are given by 

In the case of the simplest ansatz for the quark propagator, as in Eq. f)14|) . the 
Schwinger-Dyson equation (|3*4*j) reduces to an equation for the gap parameter 
A. For this reason, it is also called the gap equation. After neglecting the 
dependence of the gap on the three-momentum and after performing the 
momentum integration, one arrives at the following approximate form of the 
gap equation in Euclidean space {k^ = iko) |2"H1 !5Uj : 

A(fc 4 )^-/ ^ 4A(P4) ln ,, A , , (36) 
9tt J J p l + A 2 \h-p4\ 



where A = 2(47r) 3 / 2 fj,a~ 5 ^ 2 . If the quark wave function renormalization from 
the diagonal part of the self-energy were taken into account, the extra factor 
exp (— ^y-J would appear in the expression for A [23] • The appearance 
of the logarithm in the integrand of the gap equation (|3*B|) is an artifact of 
the long-range force in QCD, mediated by the magnetic gluon modes. An 
approximate solution to the gap equation reads 

/ 3tt 3 / 2 \ 

A(0)^Aexp(-^=). (37) 



'a. 
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The result is plotted in Fig. El Here, the value of the coupling constant 
a s = a s (n) is taken at the scale of the quark chemical potential. 
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Figure 3: The value of the gap from Eq. (|37jl as a function of the quark 
chemical potential. The coupling constant is defined at the scale of the quark 
chemical potential. The solid (dashed) line gives the result with (without) 
taking the quark wave function renormalization into account. 

It appears that, at densities that exist inside stars, i.e., 77.<10n , the 
value of the quark chemical potential cannot be much larger than 500 MeV. 
The corresponding value of the QCD coupling constant a s = a s (n) is not 
small, and, in contrast to the situation in asymptotically dense QCD, the 
Cooper pairing dynamics is not weakly interacting. This suggests that the 
above analysis of the Schwinger- Dyson equation may not be suitable for the 
description of matter inside stars. In spite of this, the solution in Eq. (|37j) is of 
fundamental importance in theory. This is one of very few non-perturbative 
solutions in QCD which is obtained from first principles, and which can be 
systematically improved by studying higher order corrections. Also, this 
solution shows that the phase diagram of QCD in the plane of temperature 
and chemical potential contains color superconducting phases at least at 
asymptotic densities. 

1.3.5 Gap equation in the NJL model 

If deconfined quark matter exists in the central regions of compact stars, it 
is likely to be color superconducting even when it is not weakly interacting. 
In order to study such a strongly interacting regime of quark matter with 
a typical value of the quark chemical potential \x ~ 500 MeV, one can use 
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various effective models of QCD. The simplest type of such models is the 
NJL model with a local four-fermion interaction in Eq. (JSJ). 

The NJL model in Eq. © lacks gluons, and, as reflection of this, it 
possesses the global instead of the local (gauge) SU(3) C color symmetry. This 
model can be viewed as a result of integrating out heavy gluons from the 
QCD action. This is possible if the gluons get nonzero masses from non- 
perturbative effects. One arrives at an effective NJL model similar to that 
in Eq. @ when the gauge fixing in QCD, needed to perform the integration, 
is consistent with the global color symmetry. 

The global color symmetry of the NJL model is broken in the 2SC phase. 
Then, because of the Goldstone theorem, five massless Nambu-Goldstone 
bosons should appear in the low-energy spectrum of such a theory. In QCD, 
however, there is no room for such Nambu-Goldstone bosons. The seeming 
contradiction is removed by noting that these Nambu-Goldstone bosons are 
not physical. Their appearance is an artifact of the gauge fixing. In particu- 
lar, there exist a gauge choice in QCD, namely the unitary gauge, in which 
these bosons can be completely eliminated. 

The gap equation in the NJL model in the mean field approximation looks 
as follows: 

A ~ ^£ [ A ( A A 

^ 'o V v /(p- /U )2 + A 2 y/(p + n) 2 + A 2 

This gap equation can be obtained from a Schwinger- Dyson equation similar 
to that in QCD in Eq. (jSlj) after the gluon long-range interaction is replaced 
by a local interaction proportional to 5{K — P). 

The approximate solution to the gap equation in Eq. reads 

A,2^^ e *p(-^ + ^). (39) 

This is very similar to the Bardeen-Cooper-Schrieffer solution in the case 
of low temperature superconductivity in solid state physics |9,. As in the 
Bardeen-Cooper-Schrieffer theory, it has the same type non-analytic depen- 
dence on the coupling constant and the same type dependence on the density 
of quasiparticle states at the Fermi surface. 

When the quark chemical potential /i takes a value in the range between 
400 MeV and 500 MeV, and the strength of the diquark pairing is Gd = r]Gs 
with r] between 0.7 and 1, the value of the gap appears to be of order 100 MeV. 
In essence, this is the result that was obtained in Ref. [T5] . 
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1.3.6 Properties of quark matter in the 2SC phase 

Now let us briefly summarize the main physical properties of the 2SC phase of 
dense quark matter. To large extent, these follow directly from the symmetry 
of the ground state, determined by the structure of the condensate in Eq. (JTUJ), 
and the magnitude of the diquark gap A. 

As we argued earlier, the blue up and blue down quarks do not participate 
in Cooper pairing in the 2SC phase. They give rise to gapless quasiparticles 
in the low energy spectrum of the theory. The density of states of such 
quasiparticles is proportional to /i 2 and, therefore, is very large. This simple 
fact has important implications for the quark matter properties. At small 
temperatures, these gapless quasiparticles give dominant contributions to 
the specific heat, as well as to the electrical and heat conductivities of the 
2SC phase. Also, the presence of the ungapped blue up and blue down 
quarks should result in a large neutrino emissivity due to the /5-processes, 
db — > u b + e~ + u e and Ub + e~ — > d b + v e FTT1. 

The other four quark quasiparticles, originating from the red and green 
paired quarks, are gapped. Their dispersion relations look like that in Eq. (jUJ) • 
At small temperatures, T < A, the contributions of these quasiparticles 
to all transport and many thermodynamic quantities are suppressed by the 
exponentially small factor exp (— A/T). Gluons do not play a very important 
role either. They are bosons and their thermal number densities are small 
at small temperatures. In addition, five out of total eight gluons are gapped 
because of the color Meissner effect. 

The presence of the unpaired blue quarks is also connected with the ab- 
sence of baryon superfluidity in the 2SC phase. The generator of the baryon 
number conservation symmetry is defined as follows: 

B = B- -|t 8 = diag color (0, 0, 1) , (40) 

where B is the generator of the U(1)b symmetry in vacuum. This vacuum 
generator B mixes with the color generator T 8 to produce the generator B 
of the U(l)s symmetry in medium. From Eq. ()40|) it is clear that only (anti- 
)blue quasiparticles carry a nonzero baryon number in the 2SC phase. 

In the 2SC phase, there is also an unbroken U(l) em gauge symmetry. The 
corresponding generator reads 

Q = Q- -±=T 8 , (41) 
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where Q = diag flavor (|, — ~) is the generator of the U(l) em symmetry in vac- 
uum. The gauge boson of U(l) em is nothing else but the medium photon. It 
should be clear, therefore, that the 2SC ground state is not subject to the 
electromagnetic Meissner effect. Thus, if there is an external magnetic field, 
it is not expelled from such a medium. 

One should notice that the quark quasiparticles carry the following Q- 
charges: 

Q(u r ) = Q{u g ) = i, Q(u b ) = l, (42) 

Q(d r ) = Q(d g ) = ~, Q(d b ) = 0. (43) 

It is interesting that the blue up quark, which remains unpaired in the 2SC 
phase, is charged (this is not true for the blue down quark which is also 
unpaired). It means that this phase of matter is a Q-conductor, and that the 
value of the electrical conductivity is very high. 

The pressure of cold dense quark matter is dominated by the Pauli pres- 
sure. The partial contribution of each quark is /i 4 /127r 2 , see Eq. (J2J). The 
correction due to color superconductivity is parametrically suppressed by 
the factor (A//i) 2 . The reason is that the Cooper pairing affects only those 
quark states that are close to the Fermi surface. In contrast, it is the whole 
Fermi sphere that contributes to the Pauli pressure. At weak coupling, the 
contribution due to diquark pairing can be easily calculated. It is given by 
(/zA/27r) 2 per each gapped quark quasiparticle in the 2SC phase, there 
are six (two flavors times three colors) quarks in total, and four of them 
give rise to quasiparticles with the gap A in their energy spectra. Thus, the 
pressure of such matter is approximately given by [321 E3] 

P^--hnr-B + 4(%-) =^~B + ^— r . 44 
127T V 2tt J 2n z tt z 

Note that the bag pressure was also added on the right hand side. If present, 
this term partially cancels the leading /i 4 term. Because of this, color super- 
conductivity may have a large effect on the hadron-quark phase transition 
[3*lj and even on the properties of compact stars [S3]. By making use of 
thermodynamic identities, one can derive the corresponding energy density, 

_ /i 2 A 2 / 2/2 8A\ 

^^6 +B+ ^[ 1 + i^)- (45) 
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The expressions in Eqs. (144)) and (|45|) give a parametric representation of the 
equation of state of dense quark matter in the 2SC phase. Such an equation 
of state is the key input in the Tolman-Oppenheimer-Volkoff equations [HE] 
which determine the interior structure of compact stars. In short, such equa- 
tions describe hydrostatic equilibrium of matter inside a star when the pres- 
sure gradient in each radial layer of the star is balanced by the gravitational 
weight of the layer itself. Several constructions of compact stars with quark 
matter in their interior were presented in Refs. |33J EH3 E3 Ell HUE HOI EE] • 

1.4 Color- flavor locked phase (Nf = 3) 

It may happen that dense baryonic matter is made not only of the lightest up 
and down quarks, but of strange quarks as well. In fact, because of a possible 
reduction in the free energy from converting non- strange quarks into strange 
quarks, one may even speculate that strange quark matter is the true ground 
state of baryonic matter 

The constituent strange quark mass in vacuum QCD is estimated to be 
of order 500 MeV. Its current mass is about 100 MeV. In dense baryonic 
matter (say, with \i ~ 500 MeV), the value of the strange quark mass should 
be in the range between 100 MeV and 500 MeV. It is plausible then that 
strange quarks also participate in Cooper pairing. 

Here, we consider an idealized version of three-flavor quark matter, in 
which all quarks are assumed to be massless. The more realistic case of a 
nonzero strange quark mass will be briefly discussed in the second lecture. 
In the massless case, the quark model possesses the global SU(3)lxSU(3).r 
chiral symmetry and the global U(l)# symmetry connected with the baryon 
number conservation. This is in addition to SU(3) C color gauge symmetry. 
Note that the generator Q = diag flavor (|, — |, — |) of the U(l) em symmetry 
of electromagnetism is traceless, and therefore it coincides with one of the 
vector-like generators of the SU(3)LxSU(3)i? chiral group. 

To large extent, the color and flavor structure of the diquark conden- 
sate of Cooper pairs in the three-flavor quark matter is again determined 
by the symmetry of the attractive diquark channel (i.e., color-antisymmetric 
antitriplet) and the Pauli exclusion principle. The spin-0 condensate corre- 
sponds to the following ground state expectation value [4"5] : 

<(^)V<>~ E c^re^ + ■ ■ ■ , (46) 
i,j=i 
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which is antisymmetric in the color and flavor indices of the constituent 
quarks, cf. Eq. (|1(J|). The 3x3 matrix Cj is determined from the global 
minimum of the free energy. It appears that Cj = 5j. The ellipsis on the 
right hand side stand for a contribution which is symmetric in color and 
flavor. A small contribution of this type is always induced in the ground 
state, despite the fact that it corresponds to a repulsive diquark channel 
|4*31 I44j . This is not surprising after noting that the symmetric condensate 
\(V' C ') 7 5 Vjf) ~ + does not break any additional symmetries |4*3*] . 

In the ground state, determined by the condensate (J46j) . the chiral sym- 
metry is broken down to its vector-like subgroup. However, the mechanism 
of this symmetry breaking is very unusual. To see this, let us rewrite the 
condensate in Eq. pUj) as follows: 

X ' 7=1 

where a, ft, a, (3 =1,2 are the spinor indices. The condensate of left-handed 
fields in Eq. (j4Tj) breaks the SU(3) C color symmetry and the SU(3)i chiral 
symmetry, but leaves the diagonal SU(3)l +c subgroup unbroken. Indeed, 
as one can check, this condensate remains invariant under the simultaneous 
flavor transformation and the compensating color transformation g c = 
(<7l) _1 . Similarly, the condensate of right-handed fields in Eq. (J47j) leaves the 
SU(3)r +c subgroup unbroken. 

When both condensates are present, the symmetry of the ground state is 
SXJ(3)l+r+ c - At the level of global symmetries, the original SU(3)lxSU(3)# 
is broken down to the vector-like SU(3)L+i?, just like in vacuum. Unlike 
in vacuum, however, this breaking does not result from any condensates 
mixing left- and right-handed fields. Instead, it results from two separate 
condensates, made of left-handed fields only and of right-handed fields only. 
The color-flavor orientations of the two condensates are "locked" to each 
other by color transformations. This mechanism is called locking, and the 
corresponding phase of matter is called color-flavor-locked (CFL) phase |4*3*] . 

To large extent, the gap equation in the three-flavor quark matter is the 
same as in the two-flavor case. The differences come only from a slightly more 
complicated color-flavor structure of the off-diagonal part of the inverse quark 
propagator (gap matrix) 



A3 = *r 



i (Ai + A 2 ) Sjfi - A 2 5tSi 



(48) 
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where two parameters Ai and A 2 determine the values of the gaps in the 
quasiparticles spectra. In the ground state that respects the SU(3)i, + ^ +c 
symmetry, the original nine quark states give rise to a singlet and an octet 
of quasiparticles. The singlet has gap Ai, and the octet has gap A 2 . 

When the small symmetric diquark condensate is neglected, one finds that 
Ai = 2A 2 , i.e., the gap of the singlet is twice as large as the gap of the octet. 
In real case, however, this relation holds only approximately. The value of 
the octet gap A 2 in asymptotically dense QCD is given by an expression 
similar to that in Eq. (|37|). but the overall coefficient is multiplied by an 
extra factor (2/3) 5/2 » 0.288 piEH]. 

1.4.1 Properties of quark matter in the CFL phase 

Let us briefly review physical properties of the CFL phase. In contrast 
to the 2SC phase, there are no gapless quark quasiparticles in the low en- 
ergy spectrum of the CFL phase. This implies that, at small temperatures, 
T< A, the contributions of quark quasiparticles to all transport and many 
thermodynamic quantities are suppressed by the exponentially small factor 
exp (— A/T). Gluons do not play any important role either. All of them are 
gapped because of the color Meissner effect. 

Unlike the 2SC phase, the CFL phase is superfluid. This is because the 
U(l)s baryon number symmetry is broken in the ground state. If such a 
phase appears in a core of a rotating star, it will be threaded with rota- 
tional vortices. This may be related to the existence of such phenomena as 
glitches. These are sudden changes of the rotational frequency observed in 
some pulsars. The glitches may be caused by occasional releasing of the an- 
gular momentum of rotational vortices, that remained pinned to the stellar 
crust l(i . 

Because of the Goldstone theorem, the breaking of the U(1)b baryon 
number symmetry should result in the appearance of a Nambu- Goldstone 
boson in the low-energy theory. In absence of gapless quark quasiparticles, 
this Nambu- Goldstone boson turns out to play an important role in many 
transport properties of cold CFL matter [33*1 147j . 

The CFL phase, like the 2SC phase, is not an electromagnetic supercon- 
ductor. Therefore, it does not expel a magnetic flux from its interior. This 
is the result of having an unbroken U(l) em gauge symmetry in the ground 
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state. The corresponding generator reads 

Q = Q - T 3 - -^T 8 . (49) 

The quark quasiparticles carry the following Q-charges: 

QM = 0, Q{u g ) = Q{u b ) = 1, (50) 
Q(d r ) = -1, Q(d g ) = Q(d b ) = 0, (51) 
Q(s r ) = -1, Q(s g ) = Q(s b ) = 0. (52) 

Since all quasiparticles are gapped, the corresponding phase of matter is a 
Q-insulator. Its electrical conductivity at small temperatures is dominated 
by thermally excited electrons and positrons |4Tj . At sufficiently low temper- 
atures, when the thermal density is low, the CFL phase becomes transparent 

to light mum. 

In order to write down the expression for the pressure of three-flavor 
quark matter in the CFL phase, we take the Pauli pressure contributions 
of nine (three flavors times three colors) quarks and add the correction due 
to color superconductivity from one quasiparticle with gap Ai and eight 
quasiparticles with gap A 2 . Thus, we derive 

p ( cfl) - ^^- fi + UrJ + H^r 

where we used the approximate relation between the singlet and octet gaps, 
Ai = 2A 2 = 2A. By making use of thermodynamic identities, one can derive 
the corresponding energy density, 

9u 4 u 2 A 2 / 2udA\ 

£(0FL) .£_ + B + 3 ^(l + ^_). (54) 

These two expressions give a parametric representation of the equation of 
state of dense quark matter in the CFL phase. 



1.4.2 Low-energy effective action 

As in the case of chiral perturbation theory in vacuum QCD one could 
write down an effective low-energy theory in the CFL phase. From the sym- 
metry breaking pattern, it is known that there are nine Nambu-Goldstone 
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bosons and one pseudo-Nambu-Goldstone boson in the low-energy spectrum 
of the theory [HUl OH]- Eight of the Nambu-Goldstone bosons are similar 
to those in vacuum QCD: three pions (it and tt^), four kaons (K°, K° 
and K^) and the eta- meson (if). The additional Nambu-Goldstone boson 
((f)) comes from breaking the U(l)s baryon symmetry. Finally, the pseudo- 
Nambu-Goldstone boson (?/) results from breaking the approximate axial 
11(1)^4 symmetry. 

The low energy action for the Nambu-Goldstone bosons in the CFL phase 
was derived in Refs. [SUl OHJ E2] • I n the chiral limit, the result reads 



C 



f 1 

J 7T 

T 
1 

+ 2 



Tr 



0^0 
'\2 2 



(d <pf-v 2 (d^ 



(55) 



where f 2 = (21 — 8 ln2)(/z/67r) 2 and v 2 = 1/3 were calculated in asymptot- 
ically dense QCD in Refs. [S3 EH] • (A nonzero mass of the t]'- meson was 
neglected here [M].) By definition, S = exp (i\ A 7r A / f^j is a unitary matrix 
field which describes the octet of the Nambu-Goldstone bosons, transforming 
under the chiral SU (3) l x SU(3)r group as follows: 



E -> UlEUI 



(56) 



where (C/ L , C/ fl ) G SU(3) L x SC/(3) fl . 

The extra terms in the low-energy action should be added when quark 
masses are nonzero EH] • Since the quark masses break the chiral symme- 
try of dense QCD explicitly, these type of corrections produce non-vanishing 
masses for all Nambu-Goldstone bosons, except one. The single boson left 
massless is the Nambu-Goldstone boson related to the baryon number break- 
ing. Of course, this symmetry is not affected by the quark masses. 

Within the framework of the high density effective theory [13 EE] , it was 
shown that an additional effect of nonzero quark masses may appear when 
the strange quark mass exceeds a critical value m s ~ ml/ 3 A 2 ^ 3 . In this case, 
the CFL phase is expected to undergo a phase transition to a phase with 
a meson (e.g., kaon or/and eta) condensate [S3 EE EH]- In the low-energy 
action of the corresponding phase, additional Nambu-Goldstone bosons with 
very unusual properties can appear 



25 



1.5 Spin-1 color superconductivity (Nf = 1) 

In the case of neutral matter in j3 equilibrium, as we shall discuss in the 
second lecture in more detail, it may happen that different quark flavors 
cannot create Cooper pairs because of a large mismatch between their Fermi 
momenta. In this case, one could consider the possibility of a much weaker 
spin-1 Cooper pairing [T2lEniEIlin21inSlEll- It is the Pauli principle that does 
not allow to construct spin-0 Cooper pairs from quarks of the same flavor. 
The color antisymmetric wave function of a pair can only be symmetric in 
spin indices. This corresponds to a spin-1 state. 

The pairing in a spin-1 triplet state and in a color antisymmetric an- 
titriplet state can lead to a rather complicated structure of the diquark con- 
densate. The general structure of the gap matrix can be written in the 
following form |64) : 

3 

A ab = iA ]T e abc C ci [k* cos 9 + T i sin o] , (57) 

c,i=l 

where k = k//c, Yj_ = Y — A; 1 (7 • k), and the explicit form of the 3x3 matrix 
C C i as well as the angular parameter 9 are determined by the minimization of 
the quark matter free energy. In the special cases with 9 = and 9 = tt/2, 
the corresponding gaps are called longitudinal and transverse, respectively. 

By choosing various specific matrices C, one could construct many differ- 
ent spin-1 phases. The four most popular of them are the so-called A-phase, 
the color-spin-locked (CSL) phase, the polar and the planar phases. These 
are determined by the following matrices 



C (A-phase) = _L , C (CSL) = 4= 1 , (58) 

^ 2 V 1 i / ^ 3 Vooiy 

/o o\ 1 / 1 \ 

C (polar) = 0, C (planar) _ I g j q ^ 



The corresponding four phases are characterized by different symmetries of 
their ground state. In particular, the original group SU(3) c xSO(3) jxU(l) em 
of one-flavor quark matter breaks down as follows [U21 EH El] : 

A-phase: SU(2) C x SO{2)j x C/(l) em) (60) 



26 



Planar 



Polar 



CSL 



SO(3)j, 

SU{2) C x SO(2)j x U(l) 
SO(2)j x U(l) em . 



(61) 
(62) 
(63) 



It was argued recently that the phase with the lowest free energy is the 
transverse A-phase 

The value of the spin-1 gap is estimated to be about two or three orders 
of magnitude smaller than a typical spin-0 gap. At best, it can be about 
1 MeV, but more realistically it is expected to be about 0.1 MeV, or less. 
The presence of such a small gap is unlikely to modify considerably any 
thermodynamic, or even transport properties of dense quark matter. It is 
fair to mention, however, that the cooling history of stars might be sensitive 
enough to feel the effect of such small gaps jHSJ. It was also suggested that 
the electromagnetic properties of spin-1 phases can be of phenomenological 
importance (53]. In contrast to the 2SC and CFL phases, discussed earlier, 
spin-1 phases reveal the electromagnetic Meissner effect. As one can see 
from Eq. (JBTj) . this clearly applies to the CSL phase. It turns out that the 
electromagnetism is also subject to the Meissner effect when there appears 
a mixture of any two independent spin-1 condensates made of quarks with 
different charges (e.g., one is made of up quarks and the other is made of 
down quarks) (HS). The type-I superconductivity in such a system could have 
observable effects on the magnetic field relaxation in pulsars. 

1.6 Summary of the first lecture 

In this lecture, only the main three color superconducting phases of dense 
quark matter have been discussed. These are the simplest possible phases 
that can be realized in one-, two-, and three-flavor quark matter. It should 
be emphasized, however, that the conditions in dense matter were idealized 
by assuming that the chemical potentials of pairing quarks were equal. In 
real situation which, for example, is realized in the interior of compact stars, 
matter happens to be very different from ideal. In particular, the charge 
neutrality and the j3 equilibrium are two very important conditions that may 
modify the quark chemical potentials and, therefore, affect the properties of 
quark matter. These issues are discussed in the second lecture. 

The aim of this lecture was to introduce the general idea of color super- 
conductivity in cold dense baryonic matter. It was argued that such matter 
is expected to be deconfined and, because of the Cooper instability, it should 
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be color superconducting. To large extent, this phenomenon is the same 
as low-temperature superconductivity in condensed matter physics [9 . One 
feature in dense QCD is very remarkable. The attractive interaction be- 
tween quarks comes from the gauge boson exchange. This is in contrast to 
the Bardeen-Cooper-Schrieffer theory in solid state physics, where the gauge 
boson (photon) exchange is repulsive, while an attractive interaction comes 
from the phonon exchange. 

One of the key points of this lecture is the observation that QCD at 
asymptotic densities is a weakly interacting (although non-perturbative) the- 
ory. Moreover, the properties of its color superconducting ground state can 
be studied analytically from first principles, providing a rare example of 
an essentially solvable limit in a non-Abelian theory. By itself, this has a 
fundamental theoretical importance. Also, this result may provide valuable 
insights in the theory of strong interactions. One of the examples might be 
the idea of duality between the hadronic and quark description of QCD [SB] . 

Although the use of weakly interacting QCD for the description of dense 
quark matter is very instructive by itself, it does not seem to be very useful for 
the quantitative description of the ground state at densities existing inside 
stars. The reason is that the density of matter in stars is bounded from 
above by the condition of the hydrostatic equilibrium, n<10no- At such not 
so large densities, one cannot treat QCD as a weakly interacting theory. In 
this situation, the use of various effective models of QCD has proved to be 
very useful. 
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2 Color superconductivity in neutral matter 



2.1 Dense matter inside stars 

As we discussed in the first lecture, it is natural to expect that some color 
superconducting phases may exist in the interior of compact stars. The esti- 
mated central densities of such stars might be sufficiently large for producing 
deconfmed quark matter. Such matter should develop the Cooper instabil- 
ity and become color superconducting. It should also be noted that typical 
temperatures inside compact stars are so low that the diquark condensate, if 
produced, would not melt. 

In the first lecture, we discussed the idealized version of dense matter, 
in which the Fermi momenta of pairing quarks were assumed to be equal. 
This does not describe the real situation that should occur inside compact 
stars. The reason is that matter in the bulk of a compact star should be 
neutral (at least, on average) with respect to electric as well as color charges. 
Otherwise, the star would not be bound by gravity which is much weaker 
than electromagnetism. Matter should also remain in (3 equilibrium, i.e., all 
P processes (e.g., such as d — » u + e~ + P e , u + e~ — > d + v e , s — > u + e~ + u e , 
and u + e~ —>■ s + u e ) should go with equal rates in both directions. 

As we shall see below, after the charge neutrality and the (3 equilibrium 
in quark matter are enforced, the chemical potentials of different quarks 
would satisfy some relations that may interfere with the dynamics of Cooper 
pairing. If this happens, some color superconducting phases may become less 
favored than others. For example, it was argued in Ref . [H7J , that a mixture 
of unpaired strange quarks and the non-strange 2SC phase, made of up and 
down quarks, is less favorable than the CFL phase after the charge neutrality 
condition is enforced. A similar conclusion was also reached in Ref. [BH] . 

Assuming that the constituent medium modified mass of the strange 
quark is large (i.e., larger than the corresponding strange quark chemical 
potential), in Ref. [HH] it was shown that neutral two-flavor quark matter 
in /3-equilibrium can have another rather unusual ground state called the 
gapless two-flavor color superconductor (g2SC). While the symmetry in the 
g2SC ground state is the same as that in the conventional 2SC phase, the 
spectrum of the fermionic quasiparticles is different. 

The existence of gapless color superconducting phases was confirmed in 
Refs. [7UJ EH H2] > and generalized to nonzero temperatures in Refs. jTHJ [71] . 
In addition, it was also shown that a gapless CFL (gCFL) phase could appear 
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in neutral strange quark matter (75J EE]- At nonzero temperature, the gCFL 
phase and several other new phases (e.g., the so-called dSC and uSC phases) 
were studied in Refs. [77J [7HJ EH] ■ If the surface tension is sufficiently small, 
as suggested in Ref. jHO], the mixed phase composed of the 2SC phase and 
the normal quark phase [SH] may be more favored. 

Here, it may be appropriate to mention that a non-relativistic analogue 
of gapless superconducting phases could appear in a trapped gas of cold 
fermionic atoms jHU IH3 Ell EH EH] • (Note that an alternative ground state 
for the atomic system, similar to the quark mixed phases in Refs. jSHl HUH IEH] , 
was proposed in Ref. |57|.) 

2.2 Gapless two-flavor color superconductivity 
2.2.1 Neutrality vs. color superconductivity 

Let us discuss the effect of charge neutrality on color superconductivity. It is 
instructive to start with the case of two-flavor quark matter first. In this case, 
the effect is most prominent and, thus, it is easiest to understand. Later, the 
results for three-flavor quark matter will be also discussed. 

In order to get an impression about the importance of the charge neu- 
trality condition in a large macroscopic chunk of matter, such as a core of a 
compact star, let us estimate the corresponding Coulomb energy. A simple 
calculation leads to the following result: 

^Coulomb ~ nl# ~ M B c* ( W ^) (j^)', (64) 

where R is the radius of the quark matter core, whose charge density is 
denoted by uq. It is easy to see that this energy is not an extensive quantity: 
the value of the corresponding energy density increases with the size of the 
system as R 2 . By taking a typical value of the charge density in the ideal 
2SC phase, nq ~ 10~ 2 e/fm 3 , the energy in Eq. ()64|) becomes a factor of 10 26 
larger than the rest mass energy of the Sun ! To avoid such an incredibly 
large energy price, the charge neutrality uq = should be satisfied with a 
very high precision. 

In the case of two-flavor quark matter, one can argue that the neutrality is 
achieved if the number density of down quarks is approximately twice as large 
as number density of up quarks, « 2n u . This follows from the fact that 
the negative charge of the down quark (Q d = —1/3) is twice as small as the 
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positive charge of the up quark (Q u = 2/3). When rid pa 2n u , the total electric 
charge density is vanishing in absence of electrons, riQ ~ Qd n d + Qu n u ~ 0. 
It turns out that even a nonzero density of electrons, required by the (3 
equilibrium condition, could not change this relation much. 

The argument goes as follows. Let us consider noninteracting massless 
quarks. In (3 equilibrium, the chemical potentials of the up quark and the 
down quark, \i u and Lid, should satisfy the relation \i d = fx u + \i e where fx e is 
the chemical potential of electrons (i.e., up to a sign, the chemical potential of 
the electric charge). By assuming that fXd ~ 2 1//3 /i u , i.e., rid ~ 2n u as required 
by the neutrality in absence of electrons, one obtains the following result for 
the electron chemical potential: /i e = fid — A*u ~ The corresponding 

density of electrons is n e pa 5 • 10~ 3 n u . We see that n e <C n u which is in 
agreement with the original assumption that rid ~ 2n„ in neutral matter. 

While the approximate relation rid pa 2n u may be slightly modified in an 
interacting system, the main conclusion would remain qualitatively the same. 
The Fermi momenta of up and down quarks, whose pairing is responsible for 
color superconductivity, are generally non-equal after the neutrality and (3 
equilibrium conditions are imposed. In this case, the Cooper pairing may 
be substantially modified, or even prevented. The study of this issue is the 
topic of the rest of this lecture. 

2.2.2 Model 

In view of strongly coupled dynamics and various non-perturbative (e.g., 
instanton) effects in dense QCD at moderate densities n<10n , existing in 
central regions of compact stars, there is no real advantage in using micro- 
scopic theory of QCD. For purposes of the current presentation, it suffices 
to use an effective model such, for example, as the NJL model with the La- 
grangian density in Eq. Q, in which the common quark chemical potential fi 
is replaced by a properly chosen color-flavor matrix fx of chemical potentials 



In f3 (chemical) equilibrium, the matrix fx is determined only by a few in- 
dependent chemical potentials which can be introduced in the corresponding 
partition function as Lagrange multipliers in front of conserved charges Qi 
of the model, 



In two-flavor quark matter, we consider only three relevant conserved charges: 



inainsiini- 




(65) 
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the baryon number the electric charge uq, and the color charge n 8 . 
Then, the matrix of quark chemical potentials is given in terms of the baryon 
chemical potential (by definition, \ib = 3/i), the electron chemical potential 
(fi e ) and the color chemical potential (/x 8 ) 



2 

Au,Q/3 — {^ij — H-eQij)$af3 + -"/=A*8^y (Tk)a/3) (66) 

where Q and T 8 are the generators of U(l) em of electromagnetism and the 
U(l)g subgroup of SU(3) C color gauge group. The explicit expressions for 
chemical potentials of different quarks read 

t-^ur f^uc 



2 1 

^ ~ 3^ e + 3^*' 


(67) 


1 1 


(68) 


^ + 3^ e + 3^' 


2 2 


(69) 




1 2 


(70) 


^ + 3^ e ~ J^ 8 ' 



Note that there exist two mutually commuting color charges in the model. 
They are related to the generators T 3 and T 8 of the SU(3) C gauge group. 
Then, generally speaking, one could introduce two independent chemical 
potentials for the corresponding two color charges. Because of the invariance 
of the 2SC ground state under the transformations of the SU(2) C color gauge 
subgroup, the introduction of the second nontrivial color chemical potential 
//3 is not necessary. This additional color chemical potential, however, is 
generally needed in three-flavor quark matter. 

Here, several comments are in order regarding the introduction of color 
quark chemical potential //§. This is required for enforcing color charge neu- 
trality only in the color superconducting ground state. In the NJL model, it 
is added by hand in the partition function, and then adjusted to neutralize 
the ground state. In QCD, on the other hand, the charge neutrality is real- 
ized dynamically due to the generation of the gluon condensate (Aq) ^ [88 . 
The appearance of such a condensate is equivalent to having a nonzero value 
of the chemical potential /i 8 ~ g s (Ao). A similar relation exists between the 
chemical potential /Z3 and the gluon condensate (Aq) in three-flavor quark 
matter. It should be clear, therefore, that nonzero gluon condensates, such as 
(Aq) or (Aq), do not break explicitly any symmetries in QCD. They appear 
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spontaneously only in the color superconducting ground state, and only after 
gauge fixing. 



2.2.3 Effective potential 



In order to derive the effective potential in the NJL model © in the mean 
field approximation, it is convenient to use the standard Hubbard-Stratonovich 
method jHH]- Then, the original NJL Lagrangian density is replaced by 
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where a, 7?, and <j) a are auxiliary fields. After taking into account the equa- 
tions of motion for these fields, 



a 



7T 



-2G S {jH>) , 
-2G S (^iisTif) , 
-2G D (t^ee a ^ 
-2G D (# C 75^e> 



c 



(72) 
(73) 
(74) 
(75) 



one can check that the Lagrangian density in (|71|) is equivalent to the original 
NJL model. In the mean field approximation, the auxiliary fields are replaced 
by their vacuum expectation values, (a) = m—mS°\ (7?) = and (0 a ) = A5g, 
and the quantum fluctuations are neglected. Then, the effective potential of 
quark matter in /3-equilibrium (with massless electrons) takes the form [73] : 



n 



+ 




+ 



(m — m^ ^ 



4Gc 



Ei + 2T In f 1 + e~ El/T ) 



(76) 



where flo is a constant added to make the pressure of the vacuum vanishing 
(the bag constant can also be included on the right hand side if necessary). 
Here, while including electrons, we neglected the contribution of neutrinos. 
This properly describes the situation inside compact stars after the delep- 
tonization had occurred, i.e., after the temperature of matter dropped so low 
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that the neutrino mean free path became larger than a typical stellar size. 
In protoneutron stars this happens in less than a minute after their birth. 

The sum in the second line of Eq. ()76|) runs over all (6 quark and 6 an- 
tiquark) quasiparticles. The dispersion relations and the degeneracy factors 
of the quasiparticles read 

E± = E(p)±» ub , [xl] (77) 

E% = E(p)±fi db , [xl] (78) 

E± ± = E±(p)±5ii. [x2] (79) 

These can be obtained in a rather straightforward way from the quark prop- 
agator. Here, the following shorthand notation was introduced: 



E(p) = ^Jp 2 + m 2 , (80) 
E±(p) = J[E(p)±jL]* + A*, (81) 



Hur + Hdg Hun + Hdr He , HS / 00 n 

" s ^^ = ^2^ = /i -y + y (82) 

r Hdg Hur Hdr Hug He /qq\ 

^ = 2 ~ 2 ~~2' [ ' 

The thermodynamic potential that determines the pressure of quark matter, 
^phys = —P, is obtained from the effective potential Q in Eq. (|75jl after 
substituting the values of Hs, He, m and A which solve the color and electric 
charge neutrality conditions, i.e., 

dQ dtt 
n 8 = -— = 0, and n Q = -— = 0, 84 

OfJLs OHe 

as well as the gap equations, i.e., 

_ = „, _ = 0. (85) 

Only a solution that satisfies Eqs. (|84j) and (J85j) can correspond to a neutral 
ground state of quark matter. If there exist several such solutions, it is the 
solution which gives the highest value of the pressure that gives the ground 
state. 
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2.2.4 Three regimes in neutral matter 

In studying neutral two-flavor quark matter, it was found that there exist 
three qualitatively different dynamical regimes which differ by the strength 
of diquark coupling, as well as by the properties of the ground state jHH [73] . 

The first regime corresponds to weak diquark coupling. In this case, a 
weak Cooper pairing is completely suppressed by the mismatch between the 
Fermi momenta of the up and down quarks. In other words, the appear- 
ance of the 2SC phase is in conflict with the constraint of charge neutrality. 
The ground state of neutral matter corresponds to the normal phase in this 
regime. One should note, though, that a much weaker spin-1 pairing be- 
tween quarks of same flavor is not forbidden in such neutral matter. In fact, 
if the temperature is sufficiently low, two independent spin-1 condensates, 
one made of up quarks and the other made of down quarks, are inevitable. 

The other limiting case is the strongly coupled regime. It is clear that, if 
the value of the diquark coupling is sufficiently large, the color condensation 
could be made as strong as needed to overcome a finite mismatch between 
the Fermi surfaces of pairing quarks. In this argument, it is taken into 
account that the charge neutrality constraint is only slightly affected by the 
value of the diquark coupling. This is the case in all models that have been 
studied so far. In this regime, the ground state is in the 2SC phase, and the 
neutrality has little effect. Here, it might be appropriate to mention that the 
unpaired blue up and blue down quarks of the 2SC phase cannot produce 
any additional spin-1 condensate. The main reason is that their chemical 
potentials are badly mismatched. 

In Fig. E| the dispersion relations of quasiparticles in the 2SC phase (i.e., 
in the strongly coupled regime) and in the normal quark matter phase (i.e., 
in the weakly coupled regime) are shown by the solid and dashed lines, re- 
spectively. The dispersion relations of unpaired blue quasiparticles are not 
shown. It should be clear, however, that they are similar to those in the 
normal phase, albeit slightly shifted in the horizontal direction. 

We see from Fig. HJ that there are two types of gapped quasiparticles in 
the spectrum of the 2SC phase. In fact, each of them is a doublet with 
respect to the unbroken SU(2) C gauge group. One of the doublets has the 
gap A + SfM, and the other has the gap A — SfM, where 5/j, = fi e /2 and A are 
the solutions of Eqs. and (J85|) that describe the ground state. When 
the mismatch parameter 5 fx is vanishing, two types of quasiparticles become 
degenerate. This corresponds to the ideal 2SC phase, considered in the first 
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Figure 4: The dispersion relations of quark quasiparticles in the 2SC phase 
(solid lines) and the normal quark matter phase (dashed lines). The low- 
energy part of the dispersion relations is shown in the right panel. 



lecture. 

The actual values of the coupling constants, which define the strongly 
and weakly coupled dynamical regimes, can differ very much from one effec- 
tive model of QCD to another. Thus, it is convenient to use the value of 
the following dimensionless ratio A/Sfj, instead of the value of the diquark 
coupling constant in a quantitative definition of such regimes. In this ratio, 
only the value of the gap A is strongly sensitive to the coupling constant. 

From our classification above, it should be clear that A/8fi = corre- 
sponds to the weakly coupled regime. This should not be confused with the 
Gd = case. In fact, in any model, there should exist a critical value 
such that A/ 5fi = when Go < G^'. As we shall see below, the strongly 
coupled regime is given by A/5fi > 1. This makes sense after noting that 
the values of the two gaps in the 2SC phase, A + 5fi and A — are positive 
when A/Sfi > 1 (see Fig. HJ). 

It turns out, that there also exists an intermediate regime, in which the 
diquark coupling is neither too weak nor too strong. In this third regime, the 
ground state is given by the so-called gapless 2SC phase [HHIIZII]- This phase 
corresponds to the ratio A/5/j, in the range from to 1 (i.e., < A/5ji < 1). 
This regime will be discussed in more detail in the next subsection. 
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2.2.5 Gapless 2SC phase 



In the case of the NJL model used in Refs. [HUES], the intermediate dynam- 
ical regime with < A/5fi < 1 is realized when the strength of the diquark 
coupling is given by 0.7<??<0.8 where r\ = Gd/Gs- To be specific, below we 
use r) = 0.75. 

Let us start from a brief discussion of the neutrality conditions in Eq. (|84j). 
At zero temperature, they are satisfied approximately when /U 8 ~ and 



It should be noted that the solutions in the two regions of parameters, A > 
6fi and A < <5/i, are very different. One could check that this reflects a 
qualitatively different nature of the quasiparticle spectra in the two regions. 
While all the modes described by the dispersion relations in Eq. (fTH|) are 
gapped when A > 5fi, there appear a gapless doublet mode when A < <5/i. 

After imposing neutrality, the chemical potentials \i e and fi 8 are not free 
parameters, but (implicit) functions of A. By making use of these functions 
He = A i e(A) and fi$ = /Us(A) in Eq. (fT6|) . one can define the effective potential 
for neutral matter as a function of the gap parameter A. This potential is 
shown graphically in Fig. (solid line). 

To emphasize the role of the neutrality condition, in Fig. the result 
for the effective potential of non-neutral quark matter at a fixed value of the 
electron chemical potential fi e ~ 148 MeV is also shown (dashed line). The 
chemical potential [i e is chosen so that the electric charge density vanishes 
when the value of the gap corresponds to the maximum of the potential. This 
effective potential (labelled "/i e = 148 MeV") describes negatively charged 
matter to the left from the maximum, and positively charged matter to 
the right from the maximum. On both sides, it is unphysical because of 
a large (infinite in an infinite volume) contribution of the Coulomb energy, 
see Eq. (164)1 . that was neglected in Eq. (fTE|) . It is fair to note that the charged 
minima of the dashed line may become physically important if mixed phases 
are allowed [HOI EHJ - 

From the location of the minimum of the effective potential for neutral 
quark matter (solid line in Fig. |5J), one determines the value of the gap 
parameter in the ground state, A « 68 MeV. It appears that A < 5fi = 



-/^, for A > <5/i, 



(86) 



A 




for A < 5/j,. 



(87) 
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Figure 5: The effective potential as a function of the diquark gap A calculated 
at a fixed value of the electron chemical potential \x e ~ 148 MeV (dashed 
line), and the effective potential defined along the neutrality line (solid line). 
The results are plotted for \i = 400 MeV and Gp = r/Gs with r\ = 0.75. 

^ e /2 « 74 MeV in such a ground state, i.e., it corresponds to the intermediate 
regime in our classification. The low energy quasiparticle spectrum in this 
phase contains additional (as compared to the 2SC phase) gapless modes 
[nniIZS]- This is also seen clearly from the analytical expressions in Eqs. 1)770- 
()79)l. Graphically, the dispersion relations are shown in Fig. El These should 
be compared with the corresponding quasiparticle dispersion relations in the 
2SC phase in Fig. H 

The most remarkable property of the quasiparticle spectrum in the g2SC 
phase is that the low energy excitations (£ C (5/i - A) are very similar to 
those in the normal phase. The only difference is that the values of the 
chemical potentials of the free up and down quarks are replaced by jj^ = 
ft ± yj (5/i) 2 — A 2 . This simple observation may suggest that the low energy 
(large distance scale) properties of the g2SC phase are similar to those in 
the normal phase [73] • For example, the Debye screening mass of the gluons 
of the SU(2) C subgroup should be nonzero. The latter, in fact, should be 
proportional to the density of the gapless modes JHTij . 

m\ D ^{ 4lBL= + 0{Sfl - A). (88) 

Note that the corresponding value for the Debye screening mass in the or- 
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Figure 6: The dispersion relations of quark quasiparticles in the g2SC phase 
(solid lines) and the normal quark matter phase (dashed lines). The low- 
energy part of the dispersion relations is shown in the right panel. 

dinary 2SC phase is vanishing [01]. The overall coefficient in Eq. (J88|) can 
be fixed by matching the value of the Debye screening mass with the known 
result in the normal phase. The Meissner screening mass of the gluons of the 
unbroken SU(2) C subgroup is vanishing, m\ M = 0, as it should be |9*U] . 

The magnetic properties of the g2SC phase are nothing like those in the 
normal phase. The values of the Meissner screening masses of five gluons 
are imaginary, indicating a chromomagnetic instability in dense two-flavor 
quark matter [00] • This instability may lead to a gluon condensation, and 
possibly to breaking the rotational or the translational symmetry in the true 
ground state of quark matter. This could be a completely new state, or this 
could be one of the states proposed earlier. For example, it may be a state 
with deformed quark Fermi surfaces [92], or it may turn into the Larkin- 
Ovchinnikov-Fulde-Ferrell state [02] , which was discussed in Ref. [01] in the 
context of quark matter. To make a certain statement, further studies in this 
direction are needed. 

2.2.6 Finite temperature properties 

The finite temperature properties of neutral quark matter in the intermediate 
coupling regime are also very unusual [731 UM- The most striking are the 
results for the temperature dependence of the gap parameter for various 
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values of the diquark coupling strengths. This is demonstrated by the plot 
in Fig. [7| As one can see, the results for not very strong coupling are very 
unusual. The gap function may take a finite value at nonzero temperatures 
even though it is exactly zero at zero temperature. This possibility comes 
about only because of the neutrality condition which introduces a mismatch 
between the Fermi momenta of pairing quarks. The thermal excitations 
smear the Fermi surfaces of the up and down quarks, and this opens the 
possibility of Cooper pairing which is forbidden at zero temperature. With 
increasing the temperature, the thermal fluctuations will eventually become 
too strong, and the pairing will be destroyed. 




10 20 30 40 50 0.7 0.74 0.78 

T (MeV) T) 



Figure 7: The temperature dependence of the diquark gap in neutral 
quark matter calculated for several values of the diquark coupling strength 
7] = Go/Gs (left panel), and the ratio of the critical temperature to the 
zero temperature gap in neutral quark matter as a function of the coupling 
strength r\ (right panel). 

By looking at the results in the left panel of Fig.0 one should immediately 
realize that the value of the ratio of the critical temperature T c to the value 
of the gap at zero temperature Ao is not a universal number. In contrast to 
the Bardeen-Copper-Schrieffer theory, the ratio T c /A in the neutral phase 
depends on the diquark coupling constant. This is shown explicitly in the 
right panel in Fig.[7| As we see, the ratio T c /A can be arbitrarily large, and 
even remain strictly infinite for a range of coupling strengths. 

2.3 Gap less color-flavor locked phase 

While considering a realistic model of strange quark matter, one should take 
into account that the strange quark mass is much larger than the masses 
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of the light up and down quarks. Most likely, the actual value of strange 
quark mass in dense medium should be in the range between about 100 MeV 
and 500 MeV. This is not negligible in comparison with a typical value of 
the quark chemical potential fi ~ 500 MeV at densities existing in compact 
stars. 

In application to the Cooper pairing dynamics responsible for color su- 
perconductivity, the main effect of a non-vanishing strange quark mass is a 
reduction of the strange quark Fermi momentum, 

fc« = ^2-m2^-|j. (89) 

The magnitude of the reduction is approximately given by the value of m 2 s /2[i. 
This quantity plays the role of a mismatch parameter in three-flavor quark 
matter, which is similar to S/i = fi e /2 in two- flavor quark matter. For obvious 
reasons, this mismatch should interfere with Cooper pairing between strange 
and non-strange quarks. 

Because of a nonzero strange quark mass, the diquark condensate of the 
ideal CFL phase (}4Tj) may get distorted as follows [HBJ HS| : 

(V>2^K) = - ~ A 1 e ijl e^ + A 2 e lj2 e^ + A 3 ^ 3 + - • • , 

(90) 

where the ellipsis denote terms symmetric in color and flavor. As in the ideal 
CFL phase, they are small and not crucial for the qualitative understanding 
of strange quark matter [TS] . 

The three gap parameters Ai, A2 and A3 in Eq. ()90|) describe down- 
strange, up-strange and up-down condensates of Cooper pairs, respectively. 
A nonzero value of rr? s j2[i affects first and foremost the pairing between the 
strange and non-strange quarks, i.e., the gap parameters Ai and A 2 . Because 
of the color-flavor locking, preserved in the diquark condensate (|90|). this 
translates into a special role played by the blue quarks. 

By starting from the massless case and gradually increasing the value of 
strange quark mass, one finds that the CFL phase stays robust until a critical 
value of the control parameter vr? s j2\i ~ A is reached [75]. Here, A is the 
gap parameter in the CFL phase. (Strictly speaking, A = Ai = A2 ~ A3 in 
the CFL phase when m s 7^ 0.) Above the critical value, the charge neutrality 
cannot be accommodated in the CFL phase, and a new gapless phase appears 
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A nice feature of the CFL phase is that it stays almost automatically 
electrically neutral [§3] • The reason is that Cooper pairing in the CFL phase 
helps to enforce the equal number densities of all three quark flavors, n u = 
nd = n s . Since the sum of the charges of the up, down and strange quarks 
add up to zero, this insures that the electric charge density is vanishing, 
n Q — \ n u — \ n <i — \n s = 0. This is exactly what happens in the CFL phase 
even at nonzero, but sub-critical values of the strange quark mass. 

It is the color rather than the electric charge neutrality that plays an 
important role in three-flavor quark matter at a nonzero value of 
It has been mentioned above that, because of the color-flavor locking, the 
blue quarks have a special status in the Cooper pairing dynamics. In order 
to avoid the violation of the color neutrality, then, a nonzero compensating 
value of the color chemical potential /x 8 — — m s/2/U is needed |H7J- (Note that, 
in QCD with dynamical gluons, a nonzero value of the gluon condensate (Aq) 
will be generated instead [EE!-) 

It was shown in Ref. [73] that the CFL phase becomes gapless when 
ml/2/j, > A. In essence, the mechanism is the same as in two- flavor quark 
matter, discussed earlier in this lecture. In order to get a slightly better 
understanding, it is instructive to consider how the strange quark mass in- 
terferes with the pairing between blue down quarks and green strange quarks 
[73]. After taking into account the shift of the strange quark Fermi momen- 
tum in Eq. ([%§]) , the effective values of the chemical potentials of blue down 
and green strange quarks become 

2 m 2 
Hdb ^ M- 3/^8- M+ 7p, (91) 

off 1 ml 2m 2 . , , 

< = " + 3" 8 -2^-17' <92) 

where we took into account that /i§ — —m -s/2/i. Up to an effective shift of the 
strange quark chemical potential, these expressions for the quark chemical 
potentials are similar to those in Eqs. (pI7)) - (f7Uj) where /i e = 0. In analogy 
with the two- flavor case, see Eq. (JH3J) . we define 

Sn = = A (93) 

as a formal value of the mismatch parameter. As soon as 5fi > A, one 
should get the gapless modes in the quasiparticle spectrum. Their dispersion 
relations are given by the same expression ()79|). except that p, — > (/^#,+/^;)/2. 
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Figure 8: The quasiparticle dispersion relations for electric and color neutral 
color superconducting quark matter at T = 0, /i = 500 MeV, and m s = 
400 MeV, see Ref. [7H]. 



A typical result for the low-energy dispersion relations of all nine quasi- 
particle are shown in Fig. |H] which is taken from Ref. [THj- It should be noted 
that the two dispersion relations in panel (d) are qualitatively the same as 
those in Fig. [HI in the gapless 2SC case. 

Physical properties of the gCFL phase are very different from those of 
the CFL phase. The presence of gapless quasiparticle modes with a density 
of states proportional to fi 2 has a large effect on the thermodynamics as well 
as on the transport. In contrast to the CFL phase which is an insulator, the 
gCFL phase is a metal with a nonzero number density of electrons. Also, the 
neutrino emissivity rate from the gCFL phase should be rather high. It is 
dominated by the /3-processes: d b — > u b + e~ + z/ e and u b + e~ — > d b + v e) where 
u b and d b denote the gapless quasiparticles whose dispersion relations are 
represented by the solid lines in panels (c) and (d) in Fig. |Hl Relatively large 
contributions to the neutrino emissivity may also come from the processes 
with additional gluons: s r , 5 — > u b + e~ + v e + #4,5,6,7 and u b + e~ — > s r , 9 + u e + 
#4,5,6, 7, where #4,5,6,7 denotes one of the gluons described by the gauge fields 
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A% A*, A% or A\. 



2.4 Mixed phases with color superconducting compo- 
nents 

In addition to homogeneous (one-component) phases, one could also study 
various mixed phases of dense quark matter. The neutrality in such phases 
is imposed not locally, but globally (i.e., on average). For example, one can 
construct a phase which is made of alternating layers of different coexisting 
components. While each of the component may have a nonzero charge den- 
sity, the average charge density of the whole construction is vanishing. The 
Coulomb energy in a mixed phase is nonzero, but it does not grow with the 
volume of the system. There is also a contribution to the total energy that 
comes from surfaces separating components of the mixed phase jSHl USB EE] ■ 
If the corresponding surface tension is sufficiently large, the mixed phase 
is disfavored energetically, and therefore it cannot appear. If, on the other 
hand, the surface tension is small, the mixed phase can be the ground state 
of matter, replacing a homogeneous gapless phase. 



2.4.1 Gibbs construction 

Here is a brief introduction to the general method of constructing mixed 
phases by imposing the Gibbs conditions of equilibrium [§Elin7|. From the 
physical point of view, the Gibbs conditions enforce the mechanical as well as 
chemical equilibrium between different components in a mixed phase. This 
is achieved by requiring that the pressure of different components inside the 
mixed phase are equal, and that the chemical potentials (fj, and /x e ) are the 
same across the whole mixed phase. For example, in application to the 
mixed phase made of the normal and the 2SC quark components |HH|, these 
conditions read 

P^W e ) = P< 2SC W B ), (94) 
p = ^ NQ) =» (2SC) , (95) 
A*. = ^ e NQ) =^ SC) . (96) 

It is easy to visualize these conditions by plotting the pressure surfaces as 
functions of chemical potentials (/i and // e ) for the components of the mixed 
phase. Graphically, this is shown in Fig. The Gibbs conditions are 
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automatically satisfied along the intersection line of the pressure surfaces 
(dark solid line in Fig. [SJ). 




Figure 9: Pressure as a function of \i = hb/3 and fi e for the hadronic phase 
(at the bottom), for the 2SC phase (on the right in front) and the normal 
quark phase (on the left). The dark thick line follows the neutrality line 
in the hadronic phase, and in two mixed phases: (i) the mixed phase of 
hadronic and normal quark matter; and (ii) the mixed phase of normal and 
color superconducting quark matter. 

Different components of the mixed phase occupy different volume frac- 
tions in space. To describe this quantitatively, it is convenient to introduce 
the volume fraction of each component. As an example, let us consider the 
quark mixed phase made of the normal and the 2SC component [SH]- We 
denote the volume fraction of the normal phase by X2SC = Vnq/V (in gen- 
eral, the notation Xb nieans "the volume fraction of phase A in a mixture 
with phase B"). Then, the volume fraction of the 2SC phase is given by 
Xnq = (1 — X2sc)- From the definition, it is clear that < Xb — 1- 

The average electric charge density of the mixed phase is determined 
by the charge densities of its components taken in the proportion of the 
corresponding volume fractions. Thus, 

< MP \», /i e ) = xS%nf 0) (* „.) + (1 - x^ Q c)nf SC Kfi, n e ). (97) 
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If the charge densities of the two components have opposite signs, one can 
impose the global charge neutrality condition, n e MP > = 0. Otherwise, a 
neutral mixed phase could not exist. In the case of quark matter, the charge 
density of the normal quark phase is negative, while the charge density of the 
2SC phase is positive along the line of the Gibbs construction (dark solid line 
in Fig. IHJ). Therefore, a neutral mixed phase exists. The volume fractions of 
its components are 



X2SC = (2SC) (ATQ)' ( 98 ) 



NQ _ n e 

(2SC) 
Tie — Tie 

V 2SC - i NQ _ n e / QQ x 

XNQ - 1 X2SC - —(NQ) (2SC)- ^ ) 

Tie Tie 

After the volume fractions have been determined from the condition of 
the global charge neutrality, the energy density of the corresponding mixed 
phase can also be calculated, 

^ (MP We) = X^ (JVQ) (/^e) + (1 - X N 2S Q c)e {2SC \»^e). (100) 

This is essentially all what one needs in order to construct the equation of 
state of the mixed phase. 



2.4.2 Surface tension and Coulomb forces 

In the above construction, the effects of the Coulomb forces and the surface 
tension between different components of the mixed phase were neglected. 
In reality, these might be important. In particular, the balance between the 
Coulomb forces and the surface tension determines the geometries of different 
components inside the mixed phase. 

In the mixed phase made of the normal and 2SC components, for example, 
nearly equal volume fractions of the two quark phases are likely to form 
alternating layers (slabs) of matter. The energy cost per unit volume to 
produce such layers scales as a 2 ^ 3 (n^ sc ^ — n^ e N ^) 2 ^ 3 where a is the surface 
tension [HE] ■ Therefore, the quark mixed ph favorable phase of matter 

only if the surface tension is not too large. The estimates of Ref. jHU] show 
that the value of a at the boundary between the normal and the 2SC phases 
is only about 5 to 10MeV/fm 2 . The maximum value of a allowed in the 
mixed phase, on the other hand, is of order 10 to 15 MeV/fm 2 when the value 
of the quark chemical potential is in the range between 400 and 500 MeV. 
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Therefore, one should conclude that the mixed phase is more favorable than 
the gapless 2SC phase. 

Here, it is fair to mention that the conclusion of Ref. [80J may not be 
final yet. The reason is that the thickness of the layers, which scales as 
a i/3^ n (2sc) _ n (NQ)^-2/3 pg| j s estimated to be of order 5 to 10 fm in the 
mixed phase [EDI EH] . This is comparable to the value of the Debye screening 
length in each of the two quark phases. Therefore, the effects of charge 
screening, which have been neglected so far and which are known to increase 
the surface energy [23, may still change the conclusion. 

2.5 Summary 

In recent years there was a lot of progress in understanding the mechanism 
of color superconductivity and in studying physical properties of various pos- 
sible color superconducting phases of dense quark matter. At present, one 
can say that the QCD theory of strong interactions at asymptotic densities 
is most likely to have a color superconducting ground state. Unfortunately, 
however, one cannot say with certainty what is the critical value of the baryon 
density at which color superconductivity first appears. Because of this, one 
does not know for sure if color superconducting phases could appear at high- 
est densities existing in the Universe. 

In nature, very dense baryonic matter exists in compact stars. Thus, 
it is possible that color superconductivity exists in the central regions of 
the heaviest compact stars. This conjecture has not yet been (dis-)proved, 
and it should come under scrutiny in the near future. In this connection, 
it might be appropriate to mention that there were recent observations on 
the cooling |lUUj of one neutron star and the radius |101j of another that led 
the authors to suggest the presence of quark matter. These suggestions have 
been disputed [1021 1103] . However, the current situation, where there is no 
unambiguous evidence that deconfined quarks play a role in compact stars, 
may well change as further observations are made. 

One of the main consequences of color superconductivity in dense matter 
is the appearance of an energy gap in the quasiparticle spectrum. In QCD, 
the typical value of the gap is 100 MeV. The presence of such a large gap 
is likely to affect many physical properties of matter, which in their turn 
should be reflected in the observational properties of stars. It is of great 
interest, therefore, to perform a systematic study of the properties of dense 
baryonic matter under conditions typical for the interior of stars. The goal 
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of such studies should be a clear theoretical prediction of the most favorable 
phases of matter inside stars. Also, one should work out the key physical 
properties of different phases that could have an unambiguous imprint on 
the observational data from stars. Such a theoretical basis would be crucial 
in an attempt to discover new (e.g., quark) states of matter in nature. 

The current understanding of color superconducting phases of QCD at 
large baryon chemical potential \i may also turn out to be useful for develop- 
ing alternative methods to treat QCD at nonzero \x on the lattice |1U4[ 111)5] . 
As is known, there are not so many methods to study QCD from first prin- 
ciples. One of them is the method of lattice calculations |10b'j . As of now, 
the lattice calculations have produced many results for QCD at finite tem- 
peratures and vanishing baryon chemical potential \x = |lU7j . There are 
also some results at nonzero but small values of the baryon chemical po- 
tential [108J. The lattice calculations at large nonzero /i are problematic 
because of the famous sign problem. At fi 7^ 0, the fermion determinant in 
the functional integral becomes complex, and the standard numerical tech- 
niques, which rely on a probability interpretation of the integrand, cannot 
be easily applied. Thus, new approaches and new ideas are needed. 

Building on the knowledge gained from lattice calculations at nonzero 
temperature as well as from the effective models of QCD at nonzero chemi- 
cal potential, one would like to get eventually a complete QCD phase diagram 
in the plane of temperature and quark chemical potential. The knowledge 
of such a diagram is not only of a theoretical importance. Different regions 
of the phase diagram describe physical conditions that can be realized, for 
example, during the evolution of the early Universe or in the heavy ion colli- 
sions. There is also a region at nonzero baryon chemical potentials which is 
responsible for the physics of stars. In order to understand the corresponding 
physics, one first has to know the locations of the critical lines and the types 
of the phase transitions. Some basic features of the QCD diagram have al- 
ready started to emerge [7SJ [HE 11091 IHU| • It will probably take many years 
before the complete phase diagram is understood in detail. 
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